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Abstract 

We study propagation of two lowest order Gaussian laser beams with different wavelengths in 
weak atmospheric turbulence. Using the Rytov approximation and assuming a slow detector we 
calculate the longitudinal and radial components of the scintillation index for a typical free space 
laser communication setup. We find the optimal configuration of the two laser beams with respect 
to the longitudinal scintillation index. We show that the value of the longitudinal scintillation for 
the optimal two-beam configuration is smaller by more than 50% compared with the value for a 
single lowest order Gaussian beam with the same total power. Furthermore, the radial scintillation 
for the optimal two-beam system is smaller by 35%-40% compared with the radial scintillation 
in the single beam case. Further insight into the reduction of intensity fluctuations is gained by 
analyzing the self- and cross-intensity contributions to the scintillation index. 
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INTRODUCTION 



Propagation of light through atmospheric turbulence is the subject of a rich and very 
active field of research owing to the many applications in free space laser communications, 



remote sensing, imaging systems and targeting O, S] . In these applications it is usually 
desirable to find ways to reduce the turbulence effects on the propagating optical beam. In 
recent years there has been a renewed interest in using partially coherent sources of light as 
a method for reducing the turbulence effects and improving the system performance. Many 
works concentrated on the case where the optical source is spatially partially coherent. 
First, the properties of the mutual coherence function of spatially partially coherent beams 
propagating in atmospheric turbulence were analyzed y,!^,^]- Later on it was shown that 
intensity fluctuations of spatially partially coherent Gaussian beams can decrease as the 
coherence of the sources decreases Additional computational evidence for the smaller 

sensitivity of spatially partially coherent beams to turbulence was obtained by considering 
the distance dependence of the coherence radius of the beam 0, 1^. More recently is 



was shown theoretically 



12| that the effective propagation distance 



ll| and experimentally 

measuring beam spreading in a turbulent medium relative to free-space spreading is larger for 
spatially partially coherent beams than for fully coherent ones. Other works demonstrated 
the improved performance obtained b y u sing spatially partially coherent beams in terms of 
the normalized intensity distribution 13 1, the mean squared beam width jl^. Ifsj]. power in 
the bucket and Strehl ratio , the scintillation index , the bit error rate (BER) 0, 
lis! and the average signal to noise ratio (SNR)[3]. 

Generation of spatially partially coherent beams in various setups, e.g., by focusing a laser 
beam on a rotating random phase screen, is quite straightforward. However, it is not clear 
to what extent such setups would fit in actual applications where reliability and compact- 
ness of the optical source play an important role. Furthermore, even though a theoretical 
method for optimizing a spatially partially coherent source was outlined in Ref. [iqI, the 
problem of implementing the method for realistic systems remains far from being resolved. 
Consequently, one has to look for other methods for generating the partial coherence, which 
are potentially more reliable and easier to optimize and to control. One very promising pos- 
sibility is to look for a temporally partially coherent input optical field generated by using 
multiple laser beams with different wavelengths. Indeed, in a typical situation in which the 
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response time of the detector is large compared with the inverse of the frequency difference 
between any pair of beams in the input field, rapidly oscillating contributions to the total 
intensity would average out. As a result, one can expect smaller values of high moments of 
the intensity compared with corresponding values in the single wavelength case. This would 
in turn result in smaller values for the scintillation index and BER and higher values for the 
SNR. 

We emphasize that generation of a temporally partially coherent source of light consisting 
of multiple laser beams with different wavelengths can be realized in a convenient and efficient 
manner by using an array of vertical external cavity surface lasers (VECSELs). These 
devices have the advantage of generating high power, high brightness wavelength tunable 



TEMqo beams (lowest order Gaussian beams) 



Individual VECSEL devices have been 



demonstrated with powers of 50 Watts running multi-mode 



2l| and up to 10 Watts, single 



221. Moreover, the semiconductor 



TEMqo mode and spectrally narrow (<0.1nm line width) 
multiple quantum well active mirror of these devices should allow for multi-GHz modulation 
rates for data transmission. 

Propagation of temporally partially coherent light in atmospheric turbulence was first 
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24| . Considering a single infinite plane wave this author obtained 



studied by Fante 

approximate analytic expressions for the scintillation index of the planar wave in the weak 
23[ and strong [24] turbulence regimes. Later studies extended these results to the cases 
of a spatially incoherent source |2^ and a source that is partially coherent in both space 
and time 126|. More recently, Kiasaleh studied the scintillation index for a multiwavelength 

nn 

infinite plane wave in weak atmospheric turbulence [23, ■ This author showed that use of 
a multiwavelength plane wave leads to an increase in the achievable SNR as well as in the 
upper bound on the SNR imposed by increasing the aperture size. These previous studies 
focused on infinite planar waves or on spherical waves, whereas in reality, Gaussian laser 
beams with finite initial spot size and possibly non-zero phase front radius of curvature 
are employed. Since the dynamics of the optical field strongly depends on the initial spot 
size and phase front radius of curvature it is important to take these characteristics into 
account when calculating turbulence effects. Furthermore, the fact that the optical field of 
Gaussian laser beams depends on 2A^ additional parameters allows for a greater flexibility 
in optimization and control of the input fleld against turbulence effects. Thus, the lower 
values for high moments of the intensity combined with the greater flexibility make the setup 
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based on multiple Gaussian beams with different wavelengths very advantageous and call 
for a detailed investigation of the optical field dynamics in this setup. 

In this paper we take this important task and study propagation of two lowest order 
Gaussian laser beams with different wavelengths in weak atmospheric turbulence. Using the 
Rytov approximation and assuming a slow detector we calculate the longitudinal component 
of the scintillation index for both the Kolmogorov spectrum and the Von Karman spectrum. 
For a typical setup of a free space laser communication system we show that the longitudinal 
scintillation exhibits a minimum as a function of the initial beam separation. We interpret 
this minimum as corresponding to the optimal configuration of the two-beam system, where 
optimization is with respect to the value of the longitudinal scintillation. Moreover, the 
longitudinal scintillation for the optimal two-beam configuration is smaller by more than 
50% compared with the longitudinal scintillation for a single lowest order Gaussian beam 
with the same total power. The existence of the longitudinal scintillation minimum is 
found to be independent of the turbulence spectrum. Furthermore, we calculate the total 
scintillation index for the optimal two-beam configuration using the Von Karman spectrum 
and find that the values of the total scintillation are smaller by 12%-64% compared with 
the corresponding values for a single Gaussian beam with the same total power. Further 
insight into this improvement is gained by decomposing the total scintillation into self- and 
cross-intensity contributions. Finally, we show that the radial scintillation for the optimal 
two-beam system is smaller by 35%-40% than the radial scintillation in the single beam case. 

The rest of the paper is organized as follows. In Section 2 we present the method used to 
calculate the scintillation index as well as general expressions for the scintillation index of 
a system consisting of N Gaussian laser beams with different wavelengths. In section 3 we 
analyze in detail the total scintillation index as well as its longitudinal and radial components 
for a two-beam system considering a typical setup of a free space laser communication 
system. Section 4 is reserved for discussion. In Appendix A we derive in detail some 
relations appearing in Section 1. 

CALCULATION OF THE SCINTILLATION INDEX 

We consider propagation of N lowest order Gaussian laser beams with different wave- 
lengths Xj, where j — 1, ...,N, in weak atmospheric turbulence. For simplicity and without 
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loss of generality we assume that the beams are linearly polarized and propagate along the 
z axis. We denote by the locations of the beam centers at the input plane 2; = 0. Thus, 
the magnitude of the total electric field E sX z = ^\^ given by 

N N 

E{r, 0,t) = J2 Ejivj, 0,t) = Y, Ujirj, 0) exp [-tu.t] , (1) 

where 



f/,(r,-,0) = Ao, X 

1 ikj 



exp 



+ ^ |r-d 



|2 



(2) 



In Eqs. (HI) and ^ Ej is the electric field due to the j-th beam, r is the radius vector 
in the xy plane, rj = r — dj, t is time, kj = 27r/Aj are the wave numbers and ujj = kjC 
are the angular frequencies, where c is the speed of light. In addition, Aqj are the initial 
on-axis amplitudes, Woj are the initial spot sizes and Eoj are the initial phase front radii of 
curvature. 

We assume that the intensity of the optical field is small enough so that its evolution is 
governed by the linear wave equation. Therefore, the dynamics of the fields Uj is described 
by 

V^Uj + k] [1 + ni(r, z)f Uj = 0, (3) 

where ni(r, z) stands for the turbulence induced fluctuation in the refractive index coefficient. 
The total optical field after propagating a distance z = L is 

N N 

E{r, L,t) = J2 Ej{r„ L,t) = J2 U,{rj, L) exp [-tu^t] . (4) 

Assuming weak turbulence, |ni(r, z)| <^ 1, Eq. Q can be approximated by 

V^Uj + k] [1 + 2ni(r, z)] U, = 0. (5) 

To solve Eq. © we use the Rytov perturbation method and express Uj as 

Uj{Tj, L) = Uoj{rj, L) exp [^^(rj, L; kj)] , (6) 

where Uoj is the field of the j-th beam in free space (in the absence of turbulence) and 
ifjj is a complex field describing the effects of turbulence on the j-th beam and calculated 
perturbatively up to the second order in ni(r, z). 



Following the approach derived in Refs. j3|, |29[ we define the beam parameters Qqj and 
Aoj, which characterize the unperturbed beams in terms of the initial phase front radii of 
curvature Fqj and spot sizes Woj, respectively 

^ A - 

We also use the beam parameters Qj and Aj, characterizing the unperturbed beams in terms 
of the phase front radii of curvature Fj and spot sizes Wj at distance L 

« =1 - — = ®0J A = = ^oj (o\ 

The time-dependent total intensity of the A^-beam system is 

N N N 

j(r, m) = Y, iji^j, ^) + E E ' (9) 

where /j = \Ej\'^ = {Ujl"^ is the intensity of the j-th beam. The intensity measured by the 
detector can be defined as the following time average 

Idetir, L) = (/(r, L, t))det = - [ dtl{v, L, t), (10) 

T Jo 

where r is the response time of the detector and {. . .)det denotes time average in the detector. 
A=.un,.ng a dow detecto. and different wavelength. A, ^ A„ fa , ^ we neglect the 
terms UjU^ j 7^ m, which are rapidly oscillating with time. [See also Ref. [22[, where a 
similar approximation was made for multiple planar laser beams] . Thus, the average of the 
total intensity of the field is given by 

AT 

Iaet{r,L) = {Iir,L,t))det ^ /,(r,-, L). (11) 

i=i 

We are interested in calculating the total scintillation index a] which is defined by 

where (...) (without any subscript) stands for average over different realizations of turbu- 
lence disorder. Using Eqs. (fTTj) and (fT^ we obtain 

a^(r,L) = ^ ^— ^ —^2 1- (13) 

(Ef=i(/.(r,,i^))^ 



Following the usual convention we decompose cr| into a longitudinal (r-independent) com- 
ponent 

a?/L)=aKO,L) (14) 

and a radial component 

a'^{r,L)^a]{r,L)-al{L). (15) 

Within the framework of the Rytov approximation the average intensity of the j-th laser 
beam is given by 0| 

L)) = lojirj, L) exp [Hij{rj, L)] , 



where 



and 



Ioj{rj,L) = exp 



2rf 
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(16) 



(17) 



Hij{rj, L) = 'iTi'^k^L / d^ / dK/«$„(fi;) x 
Jo Jo 



Iq {2A jVjC^K) exp 



(18) 



In Eq. (jlSp C, = l — z/L, K, is the wave number, $„(k) is the power spectral density of the 
refractive index fluctuations and Io{x) is the modified Bessel function of the first kind and 
zero order, /o(a;) = jQ{ix), where Jo{x) is the Bessel function of first kind and zero order 
The average of the second moment (/|(rj,L)) is given by 



(/|(r„L)) = (/,(r,.,L))2exp[ff2,(r,,L)] 



(19) 



where 



pi POO 

(k) exp 

Jo Jo 
'iK^i (1 - 6,0 



X 



Jo {2Ajrj^K) — cos 



(20) 



and Qj = 1 — Qj. The cross-intensity term {Ij{rj, L)Im{Tm, L)) is given by (see Appendix 
A) 



{Ij{rj,L)lrnirm,L)) = {Ij{rj,L)){Imirm,L)) X 

exp {i?2jm(rj, Tm] kj, km) + E2mjiSmi j'l km-, 

2Re [E3jm(rj,r^; 

^j; ^m)]} ) 



(21) 



where 



E2jmirj, Tm; kj, km) = 4:71 kjkmL dKK$„(K) Jo (k- I'yjrj - 7mrm|) X 

'o ^0 



exp 



2 y rCj hffi 



(22) 



p1 noo 

E2mj{rm,rj;km,kj) = An'^kjkmL dfi;fi;$„(fi;) Jq |7mrm - 7j*rj|) 

Jo Jo 



exp 



h h ■ 



(23) 



and 



/• 1 /"OO 

Air'^kjkmL d^ dKK^n{i^)Jo {l^llj^j - lmrm\) x 
Jo Jo 



exp 



2 y kj kfYi 



(24) 



In Eqs. (121121 7i = 1 - (©i + ^A^)^. 



SPECIFIC SETUPS, OPTIMAL CONFIGURATION AND REDUCTION OF THE 
TOTAL SCINTILLATION INDEX 

We consider in detail a system with two lowest order Gaussian beams whose wavelengths 
are Ai = 1.0 x 10~^m and A2 = 1.01 x 10~^m, and whose initial spot sizes are W^oi = W02 = 
1.0cm. The initial locations of the spot centers are di = —dy/2 and d2 = dy/2 so that 
the initial beam separation is d. We assume that the initial on-axis amplitudes are equal 
Aqi = A02 and that the beams are initially coUimated so that Fqi = Fq2 = 00. We use the 
Von Karman spectrum to describe the refractive index fluctuations 



0.033a 



2 exp(-KV/tfJ 



(25) 



where Km = 5.92//o, i^out = ^/Lq, Iq and Lq are the turbulence inner and outer scales, 
respectively, and is the refractive index structure parameter. We choose = 3.0 x 
]^Q-i5j^-2/3^ corresponding to weak atmospheric turbulence conditions, and Iq = 1.0mm, 
Lo=1.0m. We concentrate on the statistics at a propagation distance L = 1km, where the 
Rytov variance 



4 = 1.23Cy/'L''/' 



(26) 



is about 0.1 for both beams and one can indeed employ the Rytov perturbation method. 
The free space spot sizes of the beams at this distance are about 3.3cm each. 

Using the above relations we first calculate the longitudinal component of the scintillation 
index as defined by Eq. (jl4j) . The longitudinal component is of special importance when 
the spot size is large compared with the radius of the receiver's collecting lens since in this 
case the average signal to noise ratio is predominantly determined by the on-axis values of 
(/) and (/^). The dependence of a] i on the initial separation between the two beams d for 
the aforementioned values of the parameters is plotted in Fig. ^ (solid line). It can be seen 
that the curve aj i{d; L) possesses a minimum at an intermediate value of d, do = 2.8cm. 
Clearly this minimum corresponds to the optimal configuration of the two-beam system for 
the given physical parameters, where optimization is performed with respect to longitudinal 
scintillation. The figure also shows a comparison with the longitudinal scintillation of a single 
lowest order Gaussian beam with the same total power Pq = ^^o^o following two 

cases: (a) the single beam has the same initial spot size as that of each of the two beams 
(square); (b) the single beam has the same amplitude as that of each of the two beams 
(circle). One can see that the two-beam system gives a 53.4% reduction of the longitudinal 
scintillation in comparison with the single beam value in case (a) and a 56.9% reduction in 
comparison with the single beam result in case (b). The dashed line in Fig. ^corresponds 
to the result obtained by using the Kolmogorov spectrum 

$„(k) = 0.033C2k-"/1 (27) 

The dotted line corresponds to the approximate result obtained by using the Kolmogorov 
spectrum ()27|1 . assuming that for each beam rj < Wj, and expanding up to the second order 
with respect to rj/Wj. The comparison of these two results with the result based on the 
Von Karman spectrum shows that for the parameters considered here the existence of the 
minimum of cr] ^{d; L) is not very sensitive to the details of the spectrum model. In addition, 
the minimum value of the longitudinal scintillation and the location of the minimum do 
obtained by using the more realistic Von Karman spectrum are slightly smaller than the 
values obtained by using the Kolmogorov spectrum. 

In order to better understand the origin of the minimum of the longitudinal component 
of the scintillation index we further decompose a] i into a self-intensity contribution a] i ^ 
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and a cross-intensity contribution cr| ^ ^ in the following manner 

<; = ^L. + ^Mc-l, (28) 

where 

2 (/?(^y/2,L)) + (/|My/2,L)) 

= {lioJW' ^ ^ 

and 

. ,,,_ 2(/i(rfy/2,L)/2(-rfy/2,L)) 

= (mW ■ ^ ^ 

The d dependence of the two components a^ ;^ and crj i^^ is shown in Fig. |21 together with 
the d dependence of the total aj i- As expected, the cross-intensity component (jji^ is a 
decreasing function of d. In contrast, erf; ^ is increasing with increasing d due to the radial 
scintillation contribution of each beam. Thus, the existence of a minimum for aj i{d; L) is 
a direct consequence of this opposite behavior of the self- and cross-intensity contributions 
with increasing beam separation. 

When the spot size is comparable with the radius of the receiver's collecting lens the radial 
component of the scintillation index becomes important and one should take into account 
both longitudinal and radial contributions. In Figure El we present the total scintillation 
index of the two-beam system for the optimal configuration, i.e. for d = do = 2.8cm, in 
a 8cmx8cm domain centered about the z axis. The total scintillation index attains values 
in the range 0.011 < a] < 0.474. The minimum value of the scintillation is attained 
on the z axis since the two beams are already strongly overlapping at this distance, i.e., 
do/2 <Wi,W2. 

It is very instructive to analyze the reduction in the total scintillation index obtained 
by using the optimal two-beam system relative to a single lowest order Gaussian beam for 
cases (a) and (b) mentioned above. For this purpose we denote by a] ^ and ajj^ the total 
scintillation of a single beam in cases (a) and (b), respectively, and define the fractional 
reduction in the total scintillation relative to the single beam values in the two cases as 

aU r,L)-a]{ r,L) 
^L(r, L) 

and 



i?.(r,L) = ^^-"^^'7,^7;^'^^ (31) 



Rb{r, L) = . (32) 
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Figure|3]shows Ra{r, L) in the 8cm x 8cm domain. It can be seen that in case (a) the fractional 
reduction factor is in the range 0.124 < Ra < 0.639 which means that the reduction is larger 
than 12.4% everywhere within the 8cmx8cm domain and can be as large as 63.9%. Similar 
calculation shows that the relative reduction in the total scintillation in case (b) is even 
larger. In both cases, the largest improvement is obtained along the y axis and in the 
corners of the domain. The smallest improvement is along the x axis, where the total 
intensity is initially small. This behavior suggests that in a 4-beam system further increase 
of the value of R can be obtained by locating the two additional beams along the x axis. 

To get further insight into the results presented in Figs. El and HI we decompose the total 
scintillation index into self- and cross-intensity contributions a] ,, and cr|^ in the following 
manner 

a]{r, L) = a,%(r, L) + a,%(r, L) - 1, (33) 

where 

^/,s(r,L)_ (j(r,L))2 ^^^^ 

and 

Mr,L) = ^^-^^p . (35) 

The self- and cross-intensity contributions to the total scintillation index for the optimal 
two-beam system are shown in Figures El and IHl respectively. The self- intensity contribution 
attains values in the range 0.523 < a] ^ < 1.354 whereas the values of the cross-intensity 
contribution are in the range 0.052 < a] ^ < 0.582. The largest values of a] ^ are attained 
in the corners of the 8cmx8cm domain, where the a'j ^ contribution is negligible. It is also 
seen that the values of o"|^ are larger along the y axis and smaller along the x axis. This 
behavior is due to the contributions coming from the radial scintillation of each beam. The 
largest (t| ^ values are attained on the x axis. These values are comparable to the values of 
0"! g on the x axis which explains the relatively small values of Ra seen along the x axis in 
Fig. H 

In analyzing the behavior of the radial scintillation we first note that unlike the situation 
in the single-beam case, in the two-beam case cr^(r, L) is not radially symmetric. To enable 
comparison with the single-beam result we define the circularly averaged radial scintilla- 
tion in the two-beam case a^^ as the average over angle 9 of the radial scintillation index: 
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a'^j.{r,L) = {(T^{r, L))0. The r- dependence of for the optimal two-beam system is shown 
in Fig. [7| The figure also shows a comparison with the radial scintillation of a single lowest 
order Gaussian beam with the same total power and initial spot size. It can be seen that 
the value of the radial scintillation index for the optimal two-beam system at a given r is 
smaller by about 35%-40% than the corresponding value for the single beam. Therefore, 
optimizing the two-beam system with respect to the longitudinal scintillation also leads to 
a significant reduction in the radial scintillation. 

CONCLUSIONS 

We investigated the dynamics of two lowest order Gaussian laser beams with different 
wavelengths in weak atmospheric turbulence. Assuming a Von Karman turbulence spectrum 
and slow detector response and using the Rytov approximation we calculated the longitudinal 
and radial components of the scintillation index for a typical free space laser communication 
setup. We found that the longitudinal scintillation possesses a minimum as a function of the 
initial beam separation. This minimum corresponds to the ideal configuration of the two 
beams, where optimization is performed with respect to the longitudinal scintillation index. 
The longitudinal scintillation for the optimal two-beam configuration is smaller by more 
than 50% compared with the value for a single lowest order Gaussian beam with the same 
total power. We introduced the self- and cross-intensity contributions to the longitudinal 
scintillation index and explained the existence of the minimum in terms of the opposite 
behavior of these contributions with increasing beam separation. Similar calculations of 
the longitudinal scintillation with the Kolmogorov spectrum show that the existence of the 
minimum is not very sensitive to the form of the turbulence spectrum. 

In actual applications the radial component of the scintillation index might be as impor- 
tant as the longitudinal component. We therefore calculated the total scintillation index 
for the optimal two-beam configuration. We found that for the same typical setup con- 
sidered above the values of the total scintillation are smaller by 12%-64% compared with 
the values of the total scintillation for a single beam with the same total power. Further 
analysis showed that the reduction of the scintillation obtained by using the optimal two- 
beam configuration was largest along the y axis, where the centers of the two beams are 
initially located, and smallest along the x axis, where the total intensity is initially small. 
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This behavior is attributed to the relatively large values of the cross-intensity contributions 
to the total scintillation along the x axis. It also suggests that in a four-beam configuration 
the system's performance can be further improved by locating the two additional beams 
along the x-axis. Finally, we showed that the circularly averaged radial scintillation index 
for the optimal two-beam configuration is smaller by 35%-40% compared with the radial 
scintillation for a single beam with the same total power. 

The results presented in this paper open many promising pathways for future theoretical 
research. A natural extension of this study is to investigate a system with N > 2 Gaussian 
beams and to characterize the A^-dependence of the scintillation reduction relative to the 
single beam case. Another possible direction is to study the dependence of the scintillation 
index in the multi-beam multi- wavelength case on the convergence/divergence of the beams, 
i.e., on the values of the phase front radii of curvature. Optimization with respect to the 
radial scintillation index at a given radius is yet another interesting problem. From the 
applications point of view it would be very interesting to evaluate the average signal to 
noise ratio and bit error rate of the multi-beam multi-wavelength system. This will require 
calculation of the optical field after the receiver's collecting lens as well as calculation of 
aperture averaging effects. Such calculations can be carried out in a straightforward manner 



by employing the ABCD ray matrix theory |3|,|29|. 
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DERIVATION OF EQS. i^^-^ 

In this appendix we derive Eqs. (PT|l - P^ for the cross-intensity term {IjijCj, L)Im{rm, L)). 
The derivation is based on the second order Rytov approximation and follows the outline 
of the calculation of (/^(r, L)) for a single lowest order Gaussian laser beam presented in 
chapters 5 and 6 of Ref. hj]. Within the framework of the Rytov approximation the optical 
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field of tlie j-tli beam is given by Eq. Therefore, the intensity of the j-th beam is 



L) = lojirj, L) exp [^'^(rj, L; kj) + ip*{vj, L; kj 



(36) 



where Ioj{rj,L) is given by Eq. (fT7|) . It follows that 



{Ij{rj, L)/^(r„, L)) = Ioj{rj, L)/om(r„, L) exp ^ •*^*^(r, L) 



(37) 



where 



^5*^*^(r, L) = ^,(rj, L; A;,) + V'Krj, L; fc,) + ij^^^, L; A;„,) + ^^.(r^, L; fc^). (38) 



Assuming that "^f"^^ is a Gaussian random variable [3| we can use the relation 

(exp [vl/f:*V, ^)] ) = exp {(v&5r^(r, L)) + [{^f:^^\r, L)) - (^/(^^(r, L)f] } . (39) 
We expand \l/|*^*''(r, L) up to the second order with respect to ni(r, z) 



a 



vl/f:*)(r,L)^vl/f:?(r,L) + vI/(2(r,L) 



[V'ii(rj, L; fcj) + V^*i(rj, L; fc^) + ^/'mi(rm, L; k^) + ^mi(rm, ^m)] + 

L'l km) + V'rn,2(-''m) -^5 k^, 



(40) 



where the subscripts 1 and 2 in t/'i and il)2 denote first and second order. Since (rii(r, z)) = 



(vl/(*^*)(r,L)) ^ (vl/5.S(r,L)) = (^^.^(rj, L; fc,)) + (V^;2(rj, ^; ^.O) + 

m; ^m)) + {4'm2i^m: L] km)) 



(41) 



and 



(vl/(.r^(r,L))2^(vl/f:*)(r,L))2 = 



(tot) , 



(42) 



In addition. 



(vl/(*^*)^(r,L)) ^ (v&5*:f (r,L)) = {^U^^, L; k,)) + (^^^^(rj, L; A:,)) + 

(V^L(rm, Ajrrj)) + (^;;i^i(rm, L] km)) + E2jj{rj, Tj] kj, kj) + 

-^2mm (I'mj ^mi ^m; ^m) ~l~ -^2jm('"j) '^mi ^m) ~l~ -^2mj;'('^m) ^j'l ^m; ^j) ~l~ 

2Re [-E'3jm(l"j5 I"mi ^j; ^m)] ) 



(43) 
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where 



E2jm{rj, Tm] kj, km) = L; kj)t/j^^{r,:n, L; km)) 



(44) 



and 



Esjmijj, T^m] kj, km) — (^jl(rj) L] kj)'llJml(j'm, L] km))- 



(45) 



Denoting 



Eijirf kj) = (^j-2(rj, L; kj)) + -(^/'|i(rj, L; kj)) 



and using Eqs. we obtain 



(exp [^If ^(r, L)J ) = exp {2Ey(r,; k^) + 2Ei„(r^; + 

E'2jj(j^j J ^ji kj, kj) -\- E2mm{j^m^ km^ km) ~\~ E2jm(j^j ^ ^rnj kj, km)~\~ 
E2mj{^m) km, kj) + 2Re [E^jmijj, I"m) kj, km)]} ■ 



(46) 



(47) 



Noting that 



{Ij{rj, L)) = loji^j, L) exp [2Eij{Yj; kj) + E2jj{rj, r^-; kj, k 



see also Ref. j3| p. 130) we arrive at 



{Ij{Tj,L)Im{rm,L)) = {Ij{Yj,L)){Im{rm,L)) X 
exp {i?2jm(rj, Tm] kj, km) + E2mj{j^mi j'l km-, kj)-\- 

2Re [i?3jm(rj, r^; A;^, /cm)]} > 



(48) 



(49) 



which is Eq. flTTjl . 

Next we derive Eq. ()22|1 for £'2jm(rj, fm! ^m)- Within the framework of the Rytov 
perturbation method the first order term L; kj) is given by (see Ref. j^l p. 105) 



iljji{rj, L; kj) = ikj J j j dz/(K, 2) exp 



^7j(2)K-rj 



IK- 



'■7j{z){L 



2ki 



(50) 



where 



, , 1 + ajZ 



(51) 
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and 



2 i 



(52) 



In Eq. fl50|) z^(K, z) is the random amplitude of the refractive index fluctuations 



nilr, z 



J j di/(K, z) exp {iK- r) 



(53) 



K = { Kx, Ky) is the two-dimensional wave vector and k = [k,^ + Ky)^^'^. Substituting Eq. 
fl5()|l into Eq. ()44j) and using the correlation relation 



(i^(K, z)i/*(K', /)) = Fn{K, \z - /|)5(K - KOd'Kd^K' 



(54) 



where S(K) stands for the Dirac delta function and -F„(K, \z — z'\) is the two-dimensional 
spectral density of the refractive index fluctuations we obtain 

-L rL 



E2jm{rj, r^; kj, k^) = kjk^ J dz J JJ d^KF„(K, \z - z'\) x 

''lA^){L-z) ll{z'){L-z') 



IK 



exp <j iK- [7i(2)rj - -i*r,{z')Y^] - — 



kj 



kr, 



(55) 



We now change variables from z and z' io ^ = z — z' and r] = {z + z')/2. In addition, we 
use the fact that -F„(K, |/i|) is centered about ;U = to extend the integration over /i to ±cx3 
and to take z = z' = rj. This calculation yields 

-L 



pij /* f poo 

E2jm{rj,rm;kj,km) = kjkm / dr] d^K / d/ii^„(K, |/i|) x 

Jo J J J-oo 



IK 



kj 



(56) 



The two-dimensional spectral density is related to the power spectral density $ri(K) via 

1 

$„(K) = — y d/xF„(K,|/i|). (57) 

We assume that the turbulence is statistically homogeneous and isotropic so that $„(K) = 
$„(«;). Substituting Eq. ()57|) into Eq. (j^Uj) and performing integration over the angular 
coordinate in K-space we obtain 

-L 



exp 



Plj POO 

m\ kj, km) = ^rc'^kjkrn df] dK K^n{i^)Jo [k {'jjTj - 7^;^r„|) X 
Jo Jo 
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(58) 



Changing variables from t] to the normahzed distance ^ = 1 — rj/L we arrive at 




oo 



m 



Z \ Kj Km / 



(59) 



which is Eq. fl2^ . Equation for E2mj{jm^ r^-; fc^, kj) and Eq. for E^j^ivj, r^; fcj, fcm) 
are obtained in a similar manner. 
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LIST OF FIGURE CAPTIONS 

Fig. 1. Longitudinal component of the scintillation index aj i as a function of the initial 
beam separation d. The sohd hne is the result obtained by using the the Von Karman 
spectrum. The dashed line is the result obtained by using the Kolmogorov spectrum and 
the dotted line corresponds to the result obtained by using the Kolmogorov spectrum and 
assuming that rj/Wj <^ 1. The square/circle stand for the longitudinal scintillation of a 
single beam with the same total power and the same initial spot size/ amplitude, respectively. 
Fig. 2. Self- and cross-intensity contributions to the longitudinal scintillation index erf ; ^ 
and c"f_/_c' respectively, vs beam separation d. The dashed line represents cr| (^^(d; L) and the 
dotted line stands for crji^^{d; L). The sohd hne corresponds to cr] i(d; L) -\- 1/2. 
Fig. 3. Total scintillation index a'j{r,L) for the two-beam system with the optimal con- 
figuration do = 2.8cm at a propagation distance L = 1km. The figure shows a 8cmx8cm 
domain centered about the z axis. 

Fig. 4. Fractional reduction of the total scintillation index Ra{r, L) obtained by using the 
optimal two-beam system relative to a single Gaussian beam with the same total intensity 
and initial spot size. 

Fig. 5. Self-intensity contribution to the total scintillation index ^^ ^(r, L) for the optimal 

two-beam system. 

Fig. 6. Cross-intensity contribution to the total scintillation index cr|^^(r, L) for the optimal 
two-beam system. 

Fig. 7. Circularly averaged radial scintillation index for the optimal two-beam system cr^^ 
as a function of radius r (circles) . The squares correspond to the result obtained for a single 
beam with the same power and initial spot size. 
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FIG. 1: Longitudinal component of the scintillation index cr|^ as a function of the initial beam 
separation d. The solid line is the result obtained by using the the Von Karman spectrum. The 
dashed line is the result obtained by using the Kolmogorov spectrum and the dotted line corre- 
sponds to the result obtained by using the Kolmogorov spectrum and assuming that rj/Wj <^ 1. 
The square/circle stand for the longitudinal scintillation of a single beam with the same total power 
and the same initial spot size/amplitude, respectively. 
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FIG. 2: Self- and cross-intensity contributions to the longitudinal scintillation index aj ^ ^ and 
crf^p, respectively, vs beam separation d. The dashed line represents ajig{d;L) and the dotted 
line stands for aj ^ J^d; L). The solid line corresponds to crj ^{d; L) + 1/2. 
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FIG. 3: Total scintillation index aj{r,L) for the two-beam system with the optimal configuration 
do = 2.8cm at a propagation distance L = 1km. The figure shows a 8cm x 8cm domain centered 
about the z axis. 
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FIG. 4: Fractional reduction of the total scintillation index i?a(r, L) obtained by using the optimal 
two-beam system relative to a single Gaussian beam with the same total intensity and initial spot 
size. 
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FIG. 5: Self-intensity contribution to the total scintillation index aj g{r,L) for the optimal two- 
beam system. 
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FIG. 7: Circularly averaged radial scintillation index for the optimal two-beam system <T^^ as a 
function of radius r (circles). The squares correspond to the result obtained for a single beam with 
the same power and initial spot size. 
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